Close or connected? Distance and connectivity effects on transport in networks. 
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We develop an analytical approach which provides the dependence of the mean first-passage time 
(MFPT) for random walks on complex networks both on the target connectivity and on the source- 
target distance. Our approach puts forward two strongly different behaviors depending on the type 
- compact or non compact - of the random walk. In the case of non compact exploration, we show 
that the MFPT scales linearly with the inverse connectivity of the target, and is largely independent 
of the starting point. On the contrary, in the compact case the MFPT is controlled by the source- 
target distance, and we find that unexpectedly the target connectivity becomes irrelevant for remote 
targets. 

PACS numbers: 



I. INTRODUCTION 

Complex networks theory is nowadays a common tool 
to analyze a broad class of phenomena in social, biolog- 
ical or physical sciences [T]-[5]. An important issue in 
the field is to quantify the impact of topological prop- 
erties of a network on its transport properties. As a 
paradigm of transport process, random walks on com- 
plex networks have been intensely studied [IrE] , and the 
mean first-passage time (MFPT) [10] to a target node - 
which quantifies the time needed for a random walker to 
find a target on the network - has been widely used as 
an indicator of transport efficiency [Tl"Ml8| . 

A striking topological feature of many real- world com- 
plex networks is the wide distribution of the number k of 
links attached to a node - the connectivity - , as exem- 
plified by the now celebrated class of scale-free networks, 
such as internet [T!5] , biological networks [2U] > stock mar- 
kets [21] or urban traffic [22] , for which the connectivity is 
distributed according to a power law. The impact of con- 
nectivity on transport properties has been put forward in 
[H [2"5H2"S] , where it was found in different examples that 
transport towards a target node can be favored by a high 
connectivity of the target, and different functional forms 
of the dependence of the MFPT on the target connectiv- 
ity were proposed. On the other hand, the dependence of 
the MFPT on geometric properties, such as the volume 
of the network and the source to target distance, has 
been obtained recently in |26H29j . where it was shown 
that the starting position of the random walker plays a 
crucial role in the target search problem. In this con- 
text, quantifying the relative importance of distance and 
connectivity effects on transport properties on complex 
networks remains an important and widely unanswered 
question, which can be summarized as follows : is it faster 
for a random walker to find either a close, or a highly 
connected target? 

Here, we propose a general framework, applicable to a 
broad class of networks, which deciphers the dependence 
of the MFPT both on the target connectivity and on the 
source to target distance, and provides a global under- 



standing of recent results obtained on specific examples. 
Our approach highlights two strongly different behaviors 
depending on the so-called type - compact or non com- 
pact - of the random walk. In the case of non compact 
exploration, the MFPT is found to scale linearly with 
the inverse connectivity of the target, and to be widely 
independent of the starting point. On the contrary, in 
the compact case the MFPT is controlled by the source 
to target distance, and we find that unexpectedly the 
target connectivity is irrelevant for remote targets. This 
analytical approach, validated numerically on various ex- 
amples of networks, can be extended to other relevant 
first-passage observables such as splitting probabilities or 
occupations times [2"7j . 

II. MODEL AND NOTATIONS 

We are interested in the MFPT denoted T(r T |rs) of a 
discrete Markovian random walker to a target r<r, start- 
ing from a source point rg, and evolving in a network of 
N sites. We denote by A;(r) the connectivity (number of 
nearest neighbors) of site r, and by (k) its average over 
all sites with a flat measure. The corresponding degree 
distribution is denoted by p(k). We assume that at each 
time step n, the walker, at site r, jumps to one of the 
neighboring site with probability l/fc(r). Let P(r,n|r') 
be the propagator, i.e. the probability that the walker is 
at r after n steps, starting from r'. The stationary proba- 
bility distribution is then given by P s tat( r ) = k(r)/N(k), 
and it can be shown that detailed balance yields the fol- 
lowing symmetry relation : 

P(r, n|r')P st at(r') = P(v', n|r)P stat (r), (1) 

which will prove to be useful. 

We consider networks with only short range degree cor- 
relations, namely such that (fc(r)fc(r')) = (fc) 2 for |r — r'| 
larger than a cut-off distance R, where the average is 
taken over all pairs r, r' with |r — r' | fixed. This hypothe- 
sis is verified in particular by networks whose Pearson 
assortativity coefficient (30. is 0, such as Erdos-Renyi 
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networks. It is however less restrictive since local de- 
gree correlations can exist, and many networks actually 
comply with this assumption, as exemplified below. The 
hypothesis of short range degree correlations implies in 
particular that the degree distribution in a shell of radius 
r > R is identical to the degree distribution p(k) over the 
whole network, so that 

Ptat(r') ~ N r (r)/N (2) 

r'\|r— r'\=r 

where N T (r) is the number of sites r' such that |r— r'| = r. 
We then introduce the weighted average at distance r of 
a function / of two space variables defined by 

{/(r, r')}r' = j^- ]T /(r, r')P B tat(r'), (3) 

' r'/|r-r'|=r 

and the standard flat average 

(fM)" = jA* E /('>0- (4) 

rV ; r'/|r-r'|=r 

III. SCALING FORM OF THE PROPAGATOR 
FOR SCALE INVARIANT PROCESSES 

We focus hereafter on transport processes having scale 
invariant properties. It this case, we can assume that 
the propagator in the infinite network size limit after 
averaging over points at a distance r from the starting 
point, satisfies the standard scaling for |r — r'| > R : 

(Fo(r,Hr'))r«n-^/^n(-^-), (5) 

where the fractal dimension dt characterizes the accessi- 
ble volume V r oc r df within a sphere of radius r, and the 
walk dimension d w characterizes the distance r oc n 1 / dm 
covered by a random walker in n steps. A first central 
result of this paper is to show numerically that the de- 
pendence of the propagator on the connectivity of the 
target site can be actually made explicit and reads 

(P (r, n\r')) r , k cx kn^f^U (-^-) , (6) 

where the average is taken over sites r at a distance r 
from r' with fixed connectivity k. An argument support- 
ing the k dependence hypothesized in Eq. ^ is that 
it satisfies the symmetry relation of Eq. (JlJ. Numer- 
ical simulations on various examples of scale invariant 
networks, such as percolation clusters and (it, u)-flowers 
(see definition below) validate this assumption, as shown 
in Figs. [I] and [2j We stress that the scaling form Q 
is verified in the cases of both compact (d w > df) and 
non compact (d w < df) exploration. We believe that 
this result on its own can be important in the analysis of 
transport processes on networks. We show next that it 
enables to obtain the explicit dependence of first-passage 
properties on the connectivity of the target site. 




1 2 . i/d 3 

r/t » 



FIG. 1: Plot of the propagator P(rr,n|rs) for non compact 
exploration. Up: Supercritical 3D percolation networks (p = 
0.8) of different sizes, and for different fc(ry). r s is chosen in 
the center of the network, and t is small enough to avoid hits 
on the network's border. Black, red, green, blue, magenta and 
orange symbols stand respectively for k = 1, 2, 3, 4, 5 and 
6. Circles, triangles, diamonds and squares stand respectively 
for networks of size 20 3 , 25 3 , 30 3 and 40 3 . Down: (2,2,2)- 
flowers. Black, red and green circles stand respectively for 
k = 2, 6 and 18. 

IV. MEAN FIRST-PASSAGE TIME 

We now extend the theory developed in to compute 
the MFPT of a discrete Markovian random walker to a 
target r^, and obtain explicitly its dependence on fc(ry). 
As shown in [Idl the MFPT satisfies the following 
exact expression : 

P stat (r T )T(r T |r s ) = H{r T \r T ) - H(r T \r s ), (7) 

where ff(r|r') = £"=i(^(r, n|r') - P stat (r)) is the 
pseudo-Green function of the problem [33] . Note that 
averaging equation |7]) for rg covering the nearest neigh- 
bors of tt gives the expression of the averaged MFPT 
(T) Kac (r T ) expected from Kac formula [3~4l 133] : 

(T) Kac (r T ) = l/P stat (r T ) - 1 - N(k)/k(r T ) - 1, (8) 

which we will use below. 

Following we consider the large N limit of Eq. 
Making use of Eq. 0, we obtain 
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FIG. 2: Plot of the propagator P(rT,n\rs) for compact ex- 
ploration. Up: Critical 3D percolation networks (p = 0.2488) 
of different sizes, and for different fc(r-r). r 3 is chosen in the 
center of the network, and t is small enough to avoid hits on 
the network's border. Black, red, green, blue, magenta and 
orange symbols stand respectively for k = 1, 2, 3, 4, 5 and 
6. Circles, and triangles stand respectively for networks of 
size 40 3 and 50 3 . Down: (3, 3)-flowers (see [31]) for different 



fc(r-r). Black, red, 
k = 2, 4, 8 and 16. 



green, blue circles stand respectively for 



-Pstat(r T )T(r T |r s ) ~ G (r T |r T ) 



g^ (r sM . (9) 



Here Go is the usual infinite space Green function defined 
by Go(r|r') = J2^=i -fo( r , n \ r '), and ~ denotes equiva- 
lence for large N. It is useful to notice that this leading 
term of the MFPT still satisfies the Kac formula Q. We 
next take the weighted average of Eq. ^ over the source 
points and obtain: 

Pstat(r T )T rT (r) ~ G (r T |r T ) - ^ (G (r s |r T )) rs , 

(10) 

where we defined T rT (r) = {T(rr|r,s)} rs . Substituting 
the scaling ^ in Eq. ( 10 ) then yields the large N equiv- 
alence of the MFPT to a target site tt averaged over 
sources, which is valid for r > R: 



In this expression the constant Af~ depends on the con- 
nectivity k of the target and B is a constant independent 
of k ad r, which depends on the scaling function II. Wc 
now distinguish two regimes depending on the compact 
or non compact nature of the transport process, and fo- 
cus on the large r regime. 



A. Compact case d w > df 

In the compact case, d w > df, which corresponds to 
recurrent random walks, we obtain that the MFPT scales 
in the large r limit as 



T rT (r) ~ N(k)Br d ™- d i 



(12) 



This shows that unexpectedly the MFPT is asymptoti- 
cally independent of the connectivity of the target, while 
the dependence on the distance r is crucial. Eq. (11) is 
valid for r large enough (typically r > R). The depen- 
dence of Af- on k, which impacts on the MFPT for r small 
only, can be estimated by assuming that this expression 
still holds approximately for short distances. Following 
|36| . we take r = 1 in Eq. (11 1 and use the Kac formula 
d8b to obtain: 



1/fc « A k + B, 



(13) 



which provides the fc-dependence of A^. We next 
aim at evaluating B. We introduce the weighted av- 
erage of the MFPT over the tar get point r(r) = 



Er T ^tat(rT)T rT (r). Using Eq. (13 1, this quantity 



writes: 



-(r) ~ N (1 + B (k) ( r d ™- d f - 1)) 



(14) 



In the case of compact exploration, the continuous space 
limit can be defined (see [35]) and imposes r(r — > 0) = 0. 
This extra equation, based on the existence of a continu- 
ous limit, enables to evaluate B as B = l/(fc). Note that 
for fractal trees (d w — df — 1) we recover the exact result 
r(r) = Nr. Finally one has : 



(15) 



which fully elucidates the dependence of the MFPT on 
k and r. We recall here that this expression is originally 
derived for r large, and that the small r regime relies 
on the less controlled assumption that the scaling form 
of the propagator ^ holds for any distance r, and in 
particular that a continuous limit exists. It will however 
prove numerically to be accurate in various examples for 
all r values. 



B. Non compact case d w < df 

In the non compact (or transient) case, d w < df, we 
obtain that the MFPT scales in the large r limit as 



T rT (r)~iV(fc) (A k + Br d - 



(11) 



T rT (r) ~N(k)A 



k ■ 



(16) 
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This shows that the MFPT is independent of r for r 
large, as was already discussed in the literature [26] . 
The dependence on k is now fully contained in the con- 
stant Afc, which we now determine. Following |24| . we 
assume that the FPT distribution is proportional to 
exp(-Akt/(N{k))), with A = 0(1), and widely indepen- 
dent of r in agreement with the result obtained in Eq.( 16 ) 
for the first moment. This implies that the global MFPT, 
defined as the MFPT averaged over all source points and 
denoted by {T rT }, scales as {T rT } oc N(k)/k. Using the 
exact result derived in [23 : 



{Tr T } = 



H(r T \r T ) 

fstat(rT) 



(17) 



we obtain that H(tt\tt), and therefore asymptotically 
the infinite space Green function Go(r^|ry), is indepen- 
dent of k in the case of non compact exploration. This 
is checked numerically in Fig. [3] Identifying in Eq. (11) 
A-k = Go(yt\yt) /k, which is finite in the case of non 
compact exploration, we finally obtain: 



T rT (r)~N(k) 



Br d w -d s 



(18) 



As in the compact case this expression is valid for r large, 
and becomes hypothetical for r small. It reveals that 
in the case of non compact exploration, the MFPT is 
independent of r for r large, and scales as the inverse 
connectivity of the target. This behavior is in strong 
contrast with the case of compact exploration. 




FIG. 3: Numerical computation of ff (rrlrr) averaged over a 
network of a given size, as a function of the target connec- 
tivity k, on supercritical (p = 0.8) and critical (p — 0.2488) 
3D percolation network. The inset stands for the supercrit- 
ical percolation network, for three sizes 10 3 (black circles) 



15 (red triangles) and 20 (green diamonds). Equation (22 1 



gives H(rT\rr) = C. The main figure stands for the critical 
percolation network, also for three sizes (same symbols), and 
a fit in CkN d ^ /d f~ 1 (straight line). 



SUMMARY OF THE RESULTS AND 
DISCUSSION 




FIG. 4: Mean First Passage Time ((T rr (r)) fe ) for critical 
Erdos-Renyi networks, as a function of the source-target dis- 
tance r, for a various target connectivity k. Circles and tri- 
angles stand for simulation results, for two sizes of the net- 
work (1000 and 2000 nodes), straight lines stand for the zero- 
constant formula ((fc) = 2) of Eq. (191. 



Finally our central result can be summarized as fol- 
lows, where the case of marginal exploration (d w = df) 
has been obtained along the same line : 



T rT (r) 

N(k) 



1 1 

k + Jk) 
1 
k 
A 



Bln(r) 



- Br d ™- d f 



f - l) if d w > d f 



if d w = dt 



if d w < df 



(19) 



This expression is very general and shows the respective 
impact of distance and connectivity on the MFPT. In 
particular the MFPT is fully explicitly determined in the 
compact case. The positive constants A and B depend 
on the network in the case of non compact exploration. 
We comment that in both cases the target connectivity 
k plays an important role at short distances r. However 
for large source-target distances r, the fc-dependence is 
damped out in the compact case, while it remains im- 
portant in the non compact case. The r-dependence is 
found to be important in the compact case and largely 
irrelevant in the non compact case in agreement with 
previous results 26 . The question raised in introduction 
can therefore be answered as follows : in the non compact 
case connected targets are found the fastest almost inde- 
pendently of their distance, while in the compact case 
close targets are found the fastest almost independently 
of their connectivity. 

We can conclude that for self-similar networks with 
short range degree correlations, the main criterion that 
governs the behavior of T is the type (compact or non 
compact) of the random walk. In particular the existence 
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FIG. 5: Mean First Passage Time ((Tr T (r))k) as a function 
of the source-target distance r, for a various target connec- 
tivity k. Up : (2, 2)-fiowers (d w = df). Circles and trian- 
gles stand for simulation results, for two sizes of the network 
(generations 4 and 5), straight lines stand for the formula 
l/k + Bln(r) of eq. |l9j, with B = 0.24. Down: random 
(2, 2)-flowers (d w = 2.5 and df — 1.9). Circles stand for simu- 



lations results, straight lines for l/fc+l/(fc)(r 



of equation ( 19 \ 



((*> = 3) 



FIG. 6: Mean First Passage Time ((T rT (r))fe) as a function of 
the source-target distance r, for a various target connectivity 
k. Up: non compact Kozma network of size X — 400, a = 
1.0. The insight shows a translation along the y axis of A/k 
with A = 2.04 according to equation (201. As predicted, 
this quantity does not depend on k. Down: compact Kozma 
network of size X = 50, a = 2.5. The expected scaling is in 
r ' 5 : circles stand for simulation results, straight lin es stand 



for l/k + l/(fc)(r u - 5 - 1) ((k) =2.5) of equation (19 1 



of loops is irrelevant. Further comments are in order, (i) 
As stressed above, Eq. ( 19 ) is derived in the large r 
regime. Its applicability to the small r regime relies on 
the assumption that the scaling form of the propagator 
^ holds for all values of r, which is not always satisfied 
for real networks. In particular when degree correlations 
exist the relation B = l/(k) obtained in the compact case 
gives only a rough estimate, and the result of Eq. ( 19 ) 
is valid only for r larger than the correlation length, (ii) 
Our results can be extended to the case of non self-similar 
networks, still under the assumption that degree corre- 
lations are negligible. Following the method developed 
above, one can infer that 




T rT (r)~iV(fc) (A/k + g(r)) 



(20) 



where g does not depend on k and satifies g(r — > oo) = C 
in the transient case, and g(r — > oo) = oo in the recurrent 
case. The relative impact of connectivity and distance is 
therefore qualitatively the same as in the case of self- 
similar networks discussed above, (iii) Incidentally, our 
results straightforwardly yield the k dependence of the 
MFPT averaged over all source points (global MFPT). 



FIG. 7: Mean First Passage Time ((T rr (r))fe) as a function of 
the source-target distance r, for a various target connectivity 
k, on supercritical 3D percolation network [p — 0.8). For this 
network, d m ~ 2 and df = 3, the exploration is non compact. 
Circles stand for simulations results, straight lines for a fit by 
A/k + Br d ™~ d f, with A ~ 2.33 and B ~ 0.8. 



G 



We find in the large N limit : 



{Tr T } 



CN d ™/ d f 
CNlnN 
CN/k 



if d v 
if d v 
if d„ 



>d f 

= df 
<d f 



(21) 



which complements previous results obtained in 



This expression, along with Eq. (17), yields as a by- 



product the large N asymptotics of H(tt\tt)'- 



H(t t \y t ) 



ChN d ™/ d f- 

CkhxN 

C 



if d n 
if d n 
if cL 



>d f 

= df 
<d f 



(22) 



This fc-dependence of H(yt\yt) is checked numerically 
in Fig. [3] and directly validates the fc-dependence of the 
global MFPT . 



Networks of Kozma et al. - These networks, defined 
in [39] . are simple euclidian lattices in which long range 
links ("short-cuts") are added. A short-cut starts from 
each node with probability p, and leads to a node at a 
distance r where r is distributed according to a power 
law of index a. We consider here a ID euclidian lattice. 
Exploration is then compact for a > 2 and non-compact 
for a < 2. Again, Fig. [6] shows a very good agreement of 
numerical simulations with equation (19 1. 

Percolation clusters - We consider percolation clusters 
in the case of bond percolation in 3D cubic lattices. The 
critical probability is p c — 0.2488... and one has d w = 
3.88... and df = 2.58... at criticality. If p > p c , df = 3 
(euclidian dimension) and d w ~ 2. Fig. [7] shows a good 
agreement of numerical simulations with equation (19). 



VI. NUMERICAL SIMULATIONS 



VII. CONCLUSION 



We have checked our main result ( 19 1 on various ex- 



amples of networks, corresponding to compact or non 
compact random walks as detailed below. We stress that 
the zero constant formula obtained in the compact case 
is in good agreement with numerical simulations in all 
the examples that we have considered. 

Erdos-Renyi networks - Erdos-Renyi networks can be 
defined as a percolation cluster on a complete graph: for 
every pair of nodes (i,j), a link exists with probability 
p. The network is then defined as the largest cluster. We 
considered clusters at the percolation transition obtained 
for p = 1/N, for which the estimated df is 1,9 — 2, 0|37|. 
We computed numerically d w ~ 2, 9, which shows that 
exploration is compact. Numerical results of Fig. [4] are 
in very good agreement with the scaling (19). 

(u,v) -flowers - These networks are constructed recur- 
sively as described in [31] : at each step, every link is 
substituted by two paths of length u and v. We ex- 
tended this definition to (u, v, w) -flowers, for which a 
third path is added. For those networks, d w — df — 
-ln(l/u + 1/v + l/w)/ln(u) (if u < v < w). Fig. [5] 
shows a very good agreement of numerical simulations 
with equation ( 19 1, despite the small size of the networks. 

Random flowers - These networks are constructed re- 
cursively as described in [38) : at each step, every link 
is substituted by two paths of length u and v. df and 
d w are determined numerically for those networks; in 
our example (2, 2)-random flowers are compact networks 
(d w — df ~ 0,6). Fig. [5] shows a good agreement of 
numerical simulations with equation (19). 



To conclude, we have proposed a general theoretical 
framework which elucidates the connectivity and source- 
target distance dependence of the MFPT for random 
walks on networks. This approach leads to explicit solu- 
tions for self-similar networks and highlights two strongly 
different behaviors depending on the type - compact or 
non compact - of the random walk. In the case of non 
compact exploration, the MFPT is found to scale as the 
inverse connectivity of the target, and to be widely in- 
dependent of the source-target distance. On the con- 
trary, in the compact case the MFPT is controlled by the 
source-target distance, and we find thatunexpectedly the 
target connectivity is irrelevant for remote targets. The 
question raised in introduction can therefore be answered 
as follows : in the non compact case connected targets 
are found the fastest almost independently of their posi- 
tion, while in the compact case close targets are found the 
fastest almost independently of their connectivity. Last, 
we stress that following [27], this explicit determination 
of MFPTs can be straightforwardly generalized to obtain 
other relevant first-passage observables, such as splitting 
probabilities or occupation times. 
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